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The standard assumption is that all three neutrino mass states are either Dirac or Majorana. 
However, it was recently suggested by AUaverdi, Dutta and one of the authors (R.N.M.) that mixed, 
or bimodal, flavor neutrino scenarios are conceivable and are consistent with all known observations 
(these were called "schizophrenic" in the ADM paper). In that case each individual mass eigenstate 
can be either Dirac or Majorana, so that the flavor eigenstates are "large" admixtures of both. 
An example of this "bimodal" situation is to consider one mass state as a Dirac particle (with a 
sterile partner), while the other two are of Majorana type. Since only Majorana particles contribute 
to neutrinoless double beta decay, the usual dependence of this observable on the neutrino mass 
is modified within this scenario. We study this in detail and, in particular, generalize the idea 
for all possible bimodal combinations. Inevitably, radiative corrections will induce a pseudo-Dirac 
nature to the Dirac states at the one-loop level, and the effects of the pseudo-Dirac mass splitting 
will show up in the fiavor ratios of neutrinos from distant cosmological sources. Comparison of 
the effective mass in neutrinoless double beta decay as well as flavor ratios at neutrino telescopes, 
for different pseudo-Dirac cases and with their usual phenomenology, can distinguish the different 
bimodal possibilities. 

I. INTRODUCTION 

The observation of neutrino masses and mixings provides the first conclusive evidence for physics beyond the 
standard model (BSM) , so that an understanding of this phenomenon will open up one clear direction for the new 
BSM physics. Knowledge of the nature of the neutrino mass is crucial in order to make progress in this search. Unlike 
quarks and charged leptons, for which the form of the mass term is unambiguous, neutrinos are electrically neutral and 
therefore allow several possibilities. The two kinds of mass terms widely discussed are: (i) a Dirac type mass, which 
requires the theory to have a lepton number symmetry as well as a right-handed (sterile) neutrino degree of freedom; 
or (ii) a Majorana type mass, which necessitates the breaking of lepton number. Implementing the first possibility 
in the standard model requires a minimal set of assumptions, i.e., simply adding three right-handed neutrinos. In 
order to understand the small masses one requires that the associated Yukawa couplings are of order 10^ or less. 
The challenge then becomes one of understanding this tiny Yukawa coupling or at least connecting it to some other 
phenomenon that requires it. In the case of Majorana neutrinos, one can write effective dimension five operators of 
the form LHLH/M, where M represents the effect of higher scale physics. While no small couplings need be invoked 
in this case, one must understand the origin of the high scale M and explore what physics is associated with it. The 
most widely discussed theories of this type are the seesaw models where the higher scale could come from the 
breaking of new symmetries such as B ~ L or possibly a grand unified theory such as 5*0(10). 

An intermediate possibility that has also been discussed in the literature is the pseudo-Dirac scenario, where a 
tiny Majorana mass is added for either one or both of the two two-component neutrino states that make up the 
Dirac neutrino If one considers all three active neutrinos to be pseudo-Dirac, then current observations put very 
stringent constraints on the magnitude of the Majorana mass Q, i.e., ^ 10~^° eV, in order to have the pseudo-Dirac 
mass splitting small enough to remain undetected in solar neutrino oscillation experiments. Roughly speaking, for 
Vi and z/2, which contain a large amount of z/g, the pseudo-Dirac mass-squared difference should be smaller than 
E/L ~ 10~^^ eV^ for solar neutrinos, otherwise the associated oscillations would have been observed in solar neutrino 
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data. Detailed explanations can be found in Ref. This tiny splitting makes these neutrinos almost Dirac particles, 
hence the name pseudo-Dirac. 

In a recent paper [4], a new possibility for neutrino masses was pointed out, where some neutrino mass eigenstates 
are Dirac while the others are Majorana. This is only phenomenologically viable if one defines the Dirac or Majorana 
nature of neutrinos in terms of the mass eigenstates, rather than the flavor eigenstates. In this case, all neutrino 
flavors have large admixtures of both Dirac and Majorana type mass, and can be called "bimodal flavor neutrinos" 
(or schizophrenic neutrinos, as in Ref. One then needs to add as many sterile neutrino states to the standard 
model as there are Dirac mass eigenstates. This is different from the pseudo-Dirac case in the sense that the lepton 
number violating and conserving terms have comparable magnitude. Another interesting feature of bimodal neutrinos 
is that unlike the case of pseudo-Dirac flavor neutrinos, where there exist stringent constraints on the Majorana 
admixture, here the oscillations of solar neutrinos (as well as all other oscillation observations) remain unaffected. In 
other words, in conventional neutrino oscillation experiments the bimodal flavor case looks the same as the pure Dirac 
or pure Majorana case. 

An obvious place where the bimodal scenario leads to a different effect from both the pure Majorana and pseudo- 
Dirac possibilities is in the predictions for neutrinoless double beta decay. This was noted for a very specific model 
in Ref. In this paper we consider the most general implementation of this idea and present the predictions for 
neutrinoless double beta decay for the cases of both normal and inverted mass ordering. 

It was also pointed out in Ref. Q that since there is no symmetry guaranteeing the bimodal possibility, one-loop 
corrections can induce a tiny (< 10""'^^ eV) amount of Majorana mass to the mass eigenstate that had a tree level 
Dirac mass, effectively making this mass eigenstate pseudo-Dirac. Although this value is well within the constraints 
from solar neutrino observations, there are implications for astrophysical neutrinos. This is one of the new points 
explored in this paper. 

In Sec. |TT] we provide a brief review of models that could lead to the bimodal flavor neutrino scenario. Sees. IIIII 
and IIVI contain a discussion of the phenomenology of the bimodal model as it pertains to astrophysical neutrino flux 
ratios and neutrinoless double beta decay, respectively. The summary and conclusions are given in Sec. |Vl 



II. MODELS WITH ONE OR TWO DIRAC MASSES AND ONE-LOOP PSEUDO-DIRAC-NESS 

There are various possible gauge models in which the bimodal possibility for neutrinos can emerge naturally. In 
Ref. a model in which only a single mass eigenstate has a Dirac mass was considered. We briefly discuss the key 
ingredients of this model and also outline a different model in which there could be two mass eigenstates with Dirac 
masses. These models illustrate the point that the bimodal scenario leading to a large Dirac and Majorana admixture 
for flavor neutrino states can be realized within gauge models. 



A. One Dirac mass eigenstate 



Following the model in Ref. 4], an S3 symmetry is introduced, permuting the three families of SU{2) lepton 
doublets (Le, L^, L^) among themselves. This reducible representation of S3 can be decomposed as 3 = 1 + 2 so that 
the following linear combinations of lepton doublet fields, transforming as one and two dimensional representations 
of S3, turn out to be the mass eigenstates: 
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{Ll,L3) ^ (^-^(2ie-i,,-Lr),^(i^-i.)) ^2. (1) 

The 5*3 singlet field couples to the right-handed (rh) neutrino field (assumed to be an ^3 singlet) , which is isolated 
from the other two rh neutrinos by additional quantum numbers. This could either be a Zn symmetry or may even 
be the local B — L itself. For example, in Ref. [5] , the B — L quantum number is chosen such that iV^ has B — L — —5 
and Ne^T each have B — L = +A. The B — L breaking Higgs can be chosen to have quantum numbers such that only 
Ne^r have large Majorana masses (see Ref. Q for details). After integrating out the seesaw right-handed neutrinos 
Ne and Nr, the effective lepton Yukawa coupling and dimension five terms can be written as 

= hL^H^N^ + ^{L,K,f + -^{LsH^f +}i.c., (2) 
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with h, hi and /13 dimensionless coupling constants, and Hu the up-type Higgs doublet. After electroweak symmetry 
breaking, the neutrino sector has one Dirac neutrino corresponding to the mass eigenstate 1^2, two Majorana eigenstates 
Vi and i^3, as well as tribimaximal mixing (TBM) 

The Lagrangian in Eq. ([2]), together with the ^3 assignments in Eq. ([T]) and the group multiplication rules allow 
one to construct the symmetric 4x4 neutrino mass matrix in the flavor basis (z/g, v^,Vti N^), i.e.. 
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with mi = hiv^/MN^, = h^v^/Mf^^, and m2 = hvu, where u„ = (Hu). In order for the Majorana and Dirac 
mass matrix elements to have comparable magnitudes, the Yukawa coupling h must be of order 10^"'^^. This can 
be motivated in supersymmetric versions of such models, where the rh sneutrino drives inflation, and a small Dirac 
coupling is required to give consistent predictions (see Ref. ^] for details). This bimodal scenario is in contrast to 
the usual pseudo-Dirac [3| models, in which the Majorana masses are much smaller than the Dirac masses, or to the 
seesaw mechanism, in which case the Dirac masses are much smaller than the right-handed Majorana mass scale. 
The implication is that despite the large Majorana mass terms, all oscillation results remain unaffected, unlike the 
conventional pseudo-Dirac case. In fact, it is easy to show that neutrinos described by this mass matrix propagate in 
matter in the same way as those in the pure Majorana or pure Dirac case. The propagation equation contains the 
active part of M^M^, which has the same form for all the scenarios we are contemplating. 
The matrix diagonalizing Eq. ^ to diag(mi, m2, 7713, —771,2) is given by 
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where O3 = (0, 0, 0) and i?(7r/4) is the 3x3 unitary rotation matrix 
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A loop-induced pseudo-Dirac mass for 1/2 (see Sec. Ill CI) will lead to perturbations to the matrix in Eq. (jSj. In the 
simple case where the (4,4) entry of M^, is perturbed to em2, the mixing matrix Vi, is modified to 



V,' 



O3 



O3 1 M Oi^ i?(7r/4 + e/4) 



(6) 



In addition, the full Pontecorvo-Maki-Nakagawa-Sakata (PMNS) matrix will include rotations from the charged 
lepton sector. This can be described by writing down the most general Yukawa superpotential as 



— h^HdiL^cjee'' 



+ J^hrHdiLrCTeP" 



L^CTrT") + H.C. , 



(7) 



where (cTg, cr^, cr^) are gauge singlet superfields, ha (a = e, /i, r) are the Yukawa couplings, and we have assumed three 
extra Z„ symmetries that "glue" each charged lepton singlet (e'^, /i"^, r"^) to the corresponding CTq gauge singlet f^. 
The implication is that even though the neutrino mass matrix is diagonalized by the TBM matrix, there are small 
corrections from the charged lepton sector via the matrix Ui, which diagonalizes it. As a result, the final PMNS 
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matrix f/pMNS has a perturbed TBM form, if one assumes the charged lepton contributions to be smaU. The complete 
diagonahzation matrix takes the same form as Eq. ([6]), with C/tbm replaced by UjUTBM- 

In this model only the i>2 mass eigenstate has a Dirac mass, but this idea can easily be generalized in the sense that 
one or both of the other two neutrino mass eigenstates (I'l and/or 1/3) could be Dirac. The case in which two of the 
eigenstates are Dirac is motivated in the next subsection. 



B. Model with two Dirac mass eigenstates 



A model in which two Dirac mass eigenstates and one Majorana mass eigenstate appear naturally is a minimal 
B — L extension of the minimal supersymmetric standard model (MSSM), where B — L is broken by the vacuum 
expectation value (VEV) of the right-handed sneutrino [sj. Apart from the fact that there are three right-handed 
neutrinos (required for anomaly cancellation) and the gauge interactions associated with B — L, the model is essentially 
the same as the MSSM, i.e., two Higgs doublets that have zero B — L. This is therefore a minimal extension of the 
MSSM with local B — L. Furthermore, by choosing the gauge group to be SU{2)l x C/(1)/3^ x U{1)b-Li the model 
preserves gauge coupling unification with B — L breaking at the TeV scale, without any additional fields. Radiative 
corrections can allow the sneutrino fields to acquire a nonzero VEV, for certain ranges of parameters. As was noted 
in the first reference of Ref . [Sll , this leads to a neutrino mass matrix of the form 




(8) 



where the rows and columns correspond to (i^q,, z^^, T^), a = 1,2,3, and V is the superpartner of the linear combi- 
nation of B - L and I3R gauge boson, i.e. {gnYzR - GblVbl) / \/ g\ + g%L- Here gR^gsL are the gauge couplings 
of C/(l)/3j,, U{1)b-l, respectively, is the 3x3 Yukawa coupling matrix for v"^, and (z/^) is a column vector with 
components given by the three i>'^ VEVs. The parameter n is the supersymmetry (SUSY) breaking Majorana mass 
term of the B — L gaugino V. One could in fact redefine the right-handed neutrino and sneutrino states so that 
the linear combination of right-handed neutrinos that mixes with left-handed neutrinos in the Dirac mass in Eq. 
is same as the one for sneutrinos that picks up a VEV (we can call this P^). In that case, only one right-handed 
neutrino is kept in the 3x3 neutrino mass matrix in the equation above. The other two right-handed neutrinos 
remain coupled to the left-handed neutrinos only through the Yukawa couplings and are not present in Eq. ([5]). In 
the above mass matrix we have neglected small contributions that could arise from induced sneutrino VEVs, which 
can mix the right-handed neutrinos with Higgsinos. 

It is clear that this matrix leads to one linear combination of light neutrinos with Majorana mass given by the 
inverse seesaw formula Q, while the two other combinations only get a Dirac mass. With additional symmetries, 
e.g. 53 as in Ref. one could get the TBM pattern for light neutrinos. However, the main point for our discussion 
is that this model can naturally lead to one Majorana and two Dirac mass eigenstates. 

For illustration, the full symmetric 5x5 mass matrix in the flavor basis (for the TBM version where vi , 1/2 are 
Dirac type) would read 
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(9) 



In analogy to the case treated in the previous subsection, this matrix can have Majorana mass terms that are similar 
in magnitude to the Dirac masses. 

The considerations from this and the preceding subsection indicate that situations in which one or two neutrino 
mass eigenstates are Dirac and the others (or other) Majorana are possible and arise in simple models. We are 
therefore motivated to look for experimental implications of these scenarios. However, we wish to note that the small 
Dirac masses are obtained at the price of tuning the Yukawa couplings. Clearly it will be more desirable to have a 
theory that can predict these small values. 



5 



C. One-loop corrections to Dirac mass eigenstates 

Here we remark on the observation 0, [l^] that one-loop corrections to tree level Dirac state(s) pick up tiny 
Majorana corrections, if lepton number is not conserved. For simplicity, let us consider the case with one Dirac state; 
the discussion easily generalizes to the case of two Dirac states. Note that in the effective low energy Lagrangian 
of Eq. (2)), only a specific linear combination of the lepton doublets that are eigenstates of S3 appear. The charged 
lepton mass terms break this symmetry, resulting in mixings between different mass eigenstates in the finite wave 
function renormalization corrections that arise at the one-loop level. In the specific case of our first example, this will 

mean new terms of the form i5i2f'27''<9pZ^i, where S12 ^ ig^^/e ^ lO^''. Upon diagonalization of the kinetic terms, 
the new states become i^'i vi + 612^2, hence the Majorana mass term for vi in the new basis leads to a Majorana 
mass of magnitude (5^2™i ^2- The leading pseudo-Dirac contribution to the Dirac eigenstate is therefore of order 

mj^ is the atmospheric mass-squared difference), corresponding to an oscillation length 
of ~ 10 kilo parsecs (kpc). This implies that extra-galactic neutrinos from sources beyond 10 kpc will have half of 
their 1/2 component oscillate into sterile neutrinos, thus affecting the observed flavor ratios of extra-galactic neutrinos, 
which we discuss in the next section. 

We will no longer specify the magnitude of the mass-squared difference of the pseudo-Dirac neutrino Vi, but simply 
call it Smf, and analyze the phenomenological consequences in flavor ratios at neutrino telescopes and in neutrinoless 
double beta decay. 



III. EXTRA-GALACTIC NEUTRINO PHENOMENOLOGY 



Extra-galactic neutrinos by definition travel large distances in space and can have different energies, depen ding on 
their source. Thus these neutrinos can be a probe of standard and nonstandard neutrino properties (see Ref. for 
reviews). In most cases fl2] the neutrinos originate from pion (and kaon) decay, followed by muon decay (tt" — )■ ^~ -\-v^ 
and ^^-^-6^+1^6 + v^), giving the initial flavor flux ratios of $g : $° : $° = 1 : 2 : 0. However, in neutron sources 
the initial ratios are 1:0:0, with electron antineutrinos originating from /? decays 11311 : in muon-damped sources 
they become 0:1:0, since the muons (but not pious) lose energy before they decay [ij]. Although the latter two 
sources are presumably less common and harder to measure (there is less total neutrino flux) , they allow for interesting 
comparative studies with the usual pure pion source. 

In general, the initial flux composition may be described as [l^ 

($°:$° :$O)^(l:n:0). (10) 

Here the parameter n distinguishes the different types of neutrino sources: for neutron sources, the initial ratio of 
1 : : is represented by the limit n — ^ 0, whereas in muon-damped sources the initial ratio of : 1 : is the limit 
n — > 00. Pure pion sources have n = 2. In each case neutrino mixing will affect the final flavor flux ratios at Earth 
detectors, and these ratios will also depend on whether the neutrinos are pseudo-Dirac or bimodal. It is well known 
that for the initial ratios of 1 : 2 : 0, the final ratios turn out to be 1 : 1 : 1 (l2|, assuming fi — t symmetry (actually, 
it suffices to assume that 3?e(C/e3) = and 023 = n/A) and three standard neutrinos. Deviations from this symmetry 
limit will be discussed below. If some or all of the neutrinos are pseudo-Dirac, the detected flux ratios are modifled, 
and it is possible to study the effects of deviations in each different case. 

In the standard three- neutrino scenario (no pseudo-Dirac effects), the flavor conversion probability reads 

P^p = P{iy^ ^ lyp) =S^p-2j2^^ {Uo.jK,U*p^Up,) = ^ |C/„nC//3.P . (11) 

i>3 i 

However, it can be shown that if all neutrinos are pseudo-Dirac [l6j . 



PaP^J2\Ua^\'\Up.\'cOs'r■^j, (12) 

where Smf = {mf)'^ — (m")^ is the small mass-squared difference between the pseudo-Dirac pairs, and it is obvious 
that this reduces to Eq. (fTT|) for Smf = 0. In the spirit of the bimodal flavor neutrino cases discussed above, not all 
of the three states Vi could be pseudo-Dirac, but only one or two. If the corresponding SmfL/iE 1, in other words 
if L/E is large enough, the cosine term averages out to 1/2. The standard effects from neutrino mixing are therefore 
modified and neutrinos from very distant sources could probe the tiny pseudo-Dirac mass-squared differences |16i - 119i | . 
Recall that mass splittings of less than about 10^^"'^ eV^ have no effect on the solar neutrino flux 3]. 



/ SmfL 
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If one assumes that only one neutrino is pseudo-Dirac (say 1^2), then the corresponding term {i = 2) of the sum in 
Eq. p2)) is modified by a factor of 1/2, leading to the probability 

PaP = |C/alp|i7,31P + i|C/a2nC/;32|' + 1 1/«3 H I ' • (13) 

This can be extended to cases in which different combinations of neutrinos are pseudo-Dirac; the reduction factor of 
1/2 is applied to the relevant terms in each case. The measured neutrino flux, $c, is the sum of the product of each 
initial flux $° with the relevant flavor conversion probability, 

$0 = ^^/30$^ (14) 

SO that the presence of one or more pseudo-Dirac neutrinos will change the final detected fiux (and fiux ratios) 
compared to the standard case. Table U shows the observable $^/$e ratio as a function of 6*12 for the different 
combinations of pseudo-Dirac neutrinos, for /i — r symmetry and initial fluxes of 1 : 2 : 0. Note that if all three 
neutrinos are pseudo-Dirac the observed flux ratio is again 1:1, with an overall reduction in flux of 1/2. In several 
cases the ratio is independent of 612. We also observe that if 1^2 or 1^1,3 are pseudo-Dirac, then $^/$e is 1 : 1 only if 
sin^ 6*12 = |, i.e., for exact TBM. 

In the most general case one can expect deviations from exact jJ, ~ t symmetry, so that the relations ^13 — and 
^23 = f are not exact. Defining the deviation parameter 

e - J - ^23 , (15) 

the probability matrix P (without pseudo-Dirac effects) with elements Pap can be approximated as 

(1 — 2c^2Sl2 C^2'Sl2 + ^ C^2'5l2 ^ ^ \ 

5(l-c?25?2)-A W-CI2SI2) , (16) 
i(l-C?2S?2)+Ay 

where s^2 = sin^0i2, c^2 = cos^0i2 and the universal correction parameter is defined as [21I [22j 

A EEE i(2esin^26'i2 + 6*13 cos (5 sin 46*12) = 2esi2C?2 + ^6*13 cos 5 sin 46*12 . (17) 

Terms of order 0(6*13), C(e^)jand 0(^136) have been neglected in this approximation. In this case it can be shown 
that the flux ratio evolves as [2l|, 

(1 : 2 : 0) ^ (1-^2A) : (1- A) : (1- A). (18) 



TABLE I: The observed "l?,i/$e neutrino flux ratio at the detector for different combinations of pseudo-Dirac neutrinos, assuming 
the initial flux ratios of 1 : 2 : and exact /i — r symmetry. Numerical values are calculated from the global flt data in Ref. [20I ]. 
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The parameters e and A lie in the ranges —0.148 < e < 0.166 and —0.10 < A < 0.11 for the current 3a ranges [20| 
of the oscihation parameters.^ 

This general framework can be applied to the pseudo-Dirac scenario: if one or more neutrinos are pseudo-Dirac, the 
probability matrix in Eq. ()16|) will be modified, leading to different final flux ratios in each case. These probabilities 
can be written (to first order in ^13 and e) in terms of 9i2, O13, and the deviation parameters e, A, and F, where 

r = ^9i3 sin 26*12 cos (5 = 76I13S12C12 coS(5 . (19) 
o 4 

The parameter F is constrained to the range —0.026 < F < 0.026. For each case, the flavor conversion probabilities 
are given by^ 



• i/i pseudo-Dirac: 



• 1/2 and pseudo-Dirac: 



• 1^2 pseudo-Dirac: 



n';i = i(c?2^?2 + A)-F, 

PHI (1 + 3c?2 + 4e) - J A - F , (20) 



^ee — C12 +0*1 



= i + ^-^^aCl - 3c?2 + 4e) - ^A + F , (21) 
n'?-^^ J(c?2sL-A)-F, 

- + -4 + -4, 

8 8 ^^ 4 12 ' 



^ee — ''12 2 
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J^r,! = 4(c?2s?2 + A) + F, 

1 1 o „ , . . 3 



P;i = 2 - 8^12(1 + 3s?2 + 4e) - ^A + F , (22) 
Pi;^ = ^fe?2-A)-F, 



1 "Next-to-next-to-leading order" terms of second order in e and 613 have been discussed in [23l |24|| . 
^ The ratio Ptt is omitted, as it is not needed to calculate flux ratios. 



• vi and 1/3 pseudo-Dirac: 
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• 1/3 pseudo-Dirac: 
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• j^i and 1/2 pseudo-Dirac: 
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These expressions can be used to calculate the final flux ratios in each case. We discuss the most straightforwardly 
measurable flux ratio 4>^/$e and also display results for the ratio $e/<i>r, which is harder to measure. 

The plots in Figs. [1] and [2] show the variation in the flux ratios $^/<I>e and '^e/'^r with sin^0i2 for the different 
possible combinations of one or two pseudo-Dirac neutrinos, assuming the standard case of an initial flux ratio of 
1 : 2 : (a pure pion source). For comparison, the standard case without any pseudo-Dirac nature is also shown, for 
which the ratios can be approximated by 

5i«l-^«l-3A and « « 1 + 3A , (26) 

$e 1 + 2A 1-A ^ ' 

using Eq. (IT51) and neglecting quadratic terms. One can see from the plots in Figs. [T] and [5] that the two cases in 
which and V12 are pseudo-Dirac show very little dependence on 812 (compare with Table H]), even with deviations 
applied. The ratio $^/$e differs considerably from the standard case if either vi or both ui and 1^2 are pseudo-Dirac, 
and can be approximated by 

^M n o A 1 I 2 - - 3^2 2e(s^2 + ^2) 3A(3 - 4.^^ - H2) , 2r(l - 4^^^) 

dA|+ 2(1 + .2^)2 {l + sl2y 2(1 + S2,)2 + (l + sl,Y ' 
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In both cases the expressions are given to first order in the deviation parameters and the curly brackets correspond 
to the standard case [Eq. (|26l) ]. For the ratio $e/^'T, Fig.[2]shows that there are potentially strong effects if either 
or both V2 and are pseudo-Dirac, in which case 

1 13 

^«{1 + 3A} + - + -A, (29) 
n , OA1 , 8 + 4s^2- 8^12^12 -18^2-3^^2 , 3A(8 + 4sf, - Ggf^ - ^fz) , 32r(2 + gf,) 

^^,,3 ^il + 3AH- ^2 ^^,^^3 + ^2 + ,2^)3 + (2 + ,2^)3 ■ i^") 

As mentioned above, the initial flavor ratios of other interesting neutrino sources, such as neutron or niuon-damped 
sources, can be parameterised as in Eq. (fTO|) . Figures |3] and |4] indicate the dependence of the ratios $;i/<&e and <i>e/'&T 
on n for the different pseudo-Dirac combinations. It is evident that in certain cases the observed ratio can be much 
larger than in the standard case. Specifically, in the case of V2 being pseudo-Dirac, the ratios $^/$e and $e/*I'r can 
become large for n ^ oo and n — >■ 0, respectively. Expanding to first order in the deviation parameters, the ratios are 
given in these cases by 

«^oo, -P;m ^ J 1 - 0^2^12 l + c^2gi2 A l l + l + 12r 3A + 4r 3A + IGF 

(b"^ P"'^ ~ 1 9r2 «2 r,„4 4 f;s2 + a 2 2 + -, o 2 4 -, o 4 4 ' y"^^) 

•i>e refj. I ^'-12*12 ^'-12*12 J ^^12 OC]^2''12 -'^"''12*12 -'^"'-12''12 

n^(, Pi;^ ^ r i - 2cf,sL i-scf^^f^ ^) ^ 6A + 8f 3A + ler 

^01-2-1 „2 ,2 + ^4 ,4 ^ r + Q„2 ,2 0.4 „2 + a.4 ,4^ ' l^^^^ 



'I't^ Per C]^2'5l2 ''12'5l2 J 3C]^2'Sl2 9C]^2'*12 9C]^2'S 



In each case the terms in curly brackets again denote the flux ratios corresponding to the general case, without 
pseudo-Dirac neutrinos, in the same limit (n — > oo or n — > 0). Additionally, if both i^i and 1^2 are pseudo-Dirac, the 
ratio $^/3>e becomes large for n — > 00, 



O"'-^ ^ P"^ '' r 1 - r2 <i2 1 _|_ r2 2 > ^ - Af 1 

n^oo -r^M ^ I ^ '^12'^12 _ _|_J^_^12fl2 A I , ^ ^ a (oo\ 

((,•^1.2 ^ p"^i,2 ~ ] 2r2 s2 2r4 I 4r^ -2 „4 „4 ^ ' 

•i^e _fe^ ^"-12*12 ^''12''12 J 



and if both j/2 and z^3 are pseudo-Dirac, the ratio $e/$T- becomes large for n 0, 



rwo _ ^ / 1 - 2cf2gi2 , I^^£i2fi2 A \ , 1 _ 6A + 8F 3A + 16r 

^1^2.3 ^ pf2,3 -01^2-1 2 2 + 4 4 ( 3^2 2 q 4 2 + q 4 4 

-I er er I, ''12''12 ''12*12 J ^"'12''12 •^''12''12 •'"'12''12 



The plots in Figs. [3] and |4] could in principle be used to rule out certain cases. If, for instance, measurements of the 
neutrino flux ratios from a muon-damped source give $^/$e > 5, four of the six possibilities would be ruled out so 
that either 1/2 or both vi and 1^2 would have to be pseudo-Dirac neutrinos. A similar result applies for the case of 
$e/$T- > 7 and neutron sources, where only 1/2 or z/2 and could be pseudo-Dirac. 



IV. NEUTRINOLESS DOUBLE BETA DECAY PHENOMENOLOGY 

Another experimental test of the bimodal flavor neutrino scenario is neutrinoless double beta decay (see Ref. [2^ 
for reviews). In the general case with three Majorana neutrino mass eigenstates, the amplitude for this process is 
proportional to the effective Majorana mass 



E 

i=l 



= |Ci2C?3|mi| -I- S?2Cl3l"^2|e' 



-^?3l"^3|e"^| 



(35) 



with a and /3 the Majorana phases. Here the decay is mediated by light, active, and massive Majorana neutrinos, 
and one assumes that there are no other new physics contributions, such as heavy neutrino exchange, right-handed 
currents or the exchange of supersymmetric particles. 

The coherent sum (wee) contains 7 out of 9 parameters of the neutrino mass matrix and is the only observable 
carrying information about the Majorana phases. It is possible for (wee) to vanish in the case of normal neutrino 
mass ordering; this is equivalent to a zero in the (1,1) element of the low energy Majorana neutrino mass matrix. In 
the inverted ordering case {niee) cannot vanish, and the lower limit is given by 



{CI2 + 



Am^ 



17 meV, 



(36) 
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where Am\ is the mass-squared difference of atmospheric neutrinos. The final extraction of the decay half-Ufe 
is affected by uncertainties in the nuclear matrix elements so that an indisputable measurement of this process 
requires improved precision in both particle and nuclear physics parameters. Future experiments such as GERDA 
and SuperNEMO aim to reach a sensitivity of order 10 meV and should thus be able to rule out the inverted mass 
ordering, as long as the relevant parameter uncertainties are reduced. 

The standard picture of neutrinoless double beta decay is modified in the presence of pseudo-Dirac neutrinos. 

With three pseudo-Dirac neutrinos, the expression in Eq. (j35|) becomes proportional to X]i=i ^ei^^t because the 
approximately degenerate eigenstates of the pseudo-Dirac pair have opposite CP parities. This contribution is effec- 
tively vanishing {{rriee) ^ 10~^ eV) and can be neglected. However, if only one or two neutrino mass eigenstates are 
pseudo-Dirac (the bimodal scenario), one effectively has a combination of the standard case [Eq. ([55]) ] and the pure 
pseudo-Dirac case. Those neutrinos that are pseudo-Dirac do not contribute to (rriee), whereas the normal Majorana 
mass eigenstates contribute as in Eq. (|35t . 

Figs. [S] and [B] show the allowed ranges in {m^e) — X)"^' parameter space, for different combinations of pseudo-Dirac 
neutrinos and both normal and inverted neutrino mass ordering. The parameter space in the standard case is included 
for comparison. We have plotted the effective mass against the sum of masses ^m^, rather than the smallest mass 
itself, because the latter is, strictly speaking, not an observable. 

In each case, the contribution from the pseudo-Dirac pair was assumed to be vanishing so that 



N 

E 



(37) 



where the index j runs over the neutrinos that are not pseudo-Dirac, and iV = 1 or iV = 2. For instance, in the case 
where only V2 is pseudo-Dirac, the effective Majorana mass becomes 

(mee) = |c?2C?3ki| + sJajma |e*''| , (38) 

and there is only one phase, /3. One can see from the plots in Figs. [5] and [B] that in the cases of 1^2 and 1^2.3 pseudo-Dirac 
and inverted mass ordering, the lower limit for (mge) is increased by a factor of 2 4]. Explicitly, the lower bound for 
the inverted ordering becomes 



mee) « J Ami ^ ^ o"""" ~ 34 meV , (39) 



to be compared with the bound for the standard case in Eq. ([55t . 

Because of the fact that c^2 — •s^2 ~ Si2i the case in which is pseudo-Dirac results in (rriee) taking its minimal 
value in the inverted ordering. Another interesting case is when i^i 3 are pseudo-Dirac with normal ordering, where 

4 



the lower limit of {niee) is given by (Amg is the mass-squared difference of solar neutrinos) 



(mee) ^ 5^2 ^Am2 « 2.9 meV, (40) 

and the amplitude for double beta decay can never vanish, in contrast to the usual normal ordering case. 

The cases where both vi and V2 are pseudo-Dirac obviously lead to small values of (mee), since the only term 
contributing is s^3|m3|. In these cases the effective mass can lie outside the regions in which one expects it in the 
general case. Another interpretation of this would be that one of the nonstandard mechanisms of neutrinoless double 
beta decay destructively interferes with the usual mass mechanism. The strategy to test this would be to perform 
multi-isotope investigation, as the cancellation is not expected to be on the same level in different nuclei. However, 
the pseudo-Dirac suppression discussed here is the same for all nuclei. 

In summary, there are several cases for which there is a significant difference from the standard case of pure Majorana 
neutrinos. If long baseline oscillation experiments establish the neutrino ordering, and/or the neutrino mass scale 
is pinned down by cosmology or direct searches, neutrinoless double decay can distinguish the different cases. This 
illustrates the discriminative power of the process. 



V. CONCLUSION 



In conclusion, we have studied two different ways to test the bimodal (schizophrenic) neutrino hypothesis that one 
or two of the neutrino mass eigenstates are Dirac particles and the others Majorana. There are in total six nontrivial 
possible combinations, and we have performed a mostly phenomenological analysis of these scenarios. We noted that 
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(i) flux ratios of extra-galactic liigh energy neutrinos, and (ii) tlie effective mass for neutrinoless double beta decay 
are sensitive to the different possibilities, showing nonstandard behavior in many cases. Figs. [T] to [6] summarize our 
results. In brief, we found that flux ratios can differ significantly from their standard values and the effective mass can 
either lie only in certain regions or even completely outside of its standard parameter space. The examples given show 
that the many different experimental signatures provide good tests of whether neutrino masses have the bimodal (or 
pseudo-Dirac) character. We have also discussed simple beyond the standard model scenarios in which such bimodal 
features can arise. Evidence for bimodal nature of neutrino mass will require major changes in our thinking about 
the physics of neutrino mass. Indeed, the field of neutrino physics has provided many surprising results in the past, 
and the question of neutrino mass origin is far from settled. If the hypothesis of bimodal neutrinos is supported by 
the experiments outlined here, it will not only provide a major departure from our current thinking about the nature 
of neutrino masses but also its theoretical origin from physics beyond the standard model. As such it will have major 
impact on the physics at and beyond the TeV scale. 
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FIG. 1: The observable flux ratio <l?^/"l>e against sin^0i2, assuming the initial neutrino flux ratios of 1 : 2 : and different 
combinations of pseudo-Dirac neutrinos (denoted by red dashed lines), with the parameters ^13, 623, and 5 varying in their 
allowed So" ranges. The black hatched region shows the general case with no pseudo-Dirac neutrinos; the red cross (black plus 
sign) shows the value of <l?f,./"l>e in each pseudo-Dirac case (the general case), assuming TBM. 



13 




0.5 1 



0.28 0.3 0.32 0.34 0.36 0.28 0.3 0.32 0.34 0.36 



sm 6 



12 



FIG. 2: Same as Fig. [T] for the observable flux ratio 
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FIG. 3: The observable flux ratio $^/$e against n, assuming the initial neutrino flux ratios of 1 : n : and different combinations 
of pseudo-Dirac neutrinos (denoted by red dashed-dotted lines), with the parameters 623, and 5 varying in their allowed 
3(7 ranges. The black hatched region shows the general case with no pseudo-Dirac neutrinos; the green dashed line shows the 
value of $;j/"l>e in each pseudo-Dirac case, assuming TBM. 
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FIG. 4: Same as Fig. |3] for the observable flux ratio 
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FIG. 5: Allowed regions in the (mee) — plane for the three different cases of one pseudo-Dirac neutrino (indicated on the 

right of each row). The black regions are for exact TBM, and the light red (green) shaded regions correspond to the 3a ranges 
of the oscillation parameters for normal (inverted) ordering. The solid (dashed) lines indicate the best-fit (3a) allowed regions 
in the standard three-neutrino scenario. 
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FIG. 6: Same as Fig. [5] for the three different cases of two pseudo-Dirac neutrinos (indicated on the right of each row). 



